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Section 1.0 


m... 


SUMMARY 


The final report on the Electrical Isolation Analysis was com- 
pleted and is included as Appendix' A. Several full size 
qualification modules were constructed using polyurethane 
encapsulant and a Tedlar front cover. 
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Boction 2,0 


INTRODUCTION 




This pi‘ogram will develop the analytical methodology for advanced 
encapsulation designs. Prom these methods design sensitivities 
will be established for the development of photovoltaic module 
criteria and the definition of needed research tasks. 


The program consists of four phases. J.n Phase X analytical models 
were developed to perform optical, thermal, electrical and struc- 
tural analyses on candidate encapsulation systems. From these 
analyses several candidate encapsulation systems were selected 
for qualification testing during Phase II. Additionally, during 
Phase II test specimens of various types will be constructed and 
tested to determine the validity of the analysis methodology 
developed in Phase X. 

During Phase 111 the following items will be covered: 

1. Coxu’cction of identified deficiencies and/or discrepancies 
between analytical models developed during Phase I and rele- 
vant test data obtained during Phase li of the above con- 
tract. 

2. Improvement and extension of prediction capability of present 
analytical models. 




3. Generation of encapsulation engineering generalities, prin- 
ciples, and design aids for photovoltaic module design. 

From these items the sensitivity of module performance to various 
material properties will be determined. This study will enable 
the intelligent direction of research into assessment of module 
life potential by analyzing those materials and their properties 
which through aging would most influence module performance. 

In Phase IV a final optimum design based on knowledge gained in 
Phases I, II and III will be developed and delivered to JPL. 
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Section 3.0 


TECHNICAL DISCUSSION 


^ 1 Electrical Report 


NASTRAN is a comprehensive finite element computer program nor- 
mally used to analyze the structural and, thermal aspects of 
complex systems. , One of the attractive features of NASTRAN is 
the ability to quickly and easily construct complex 3-D models. 
Electric fields are analogous with thermal fields and can be 
calculated with a thermal analyzer program such as NASTRAN. 
Several analyses have been done using this technique. The final 
report on the electrical analyses is presented as Appendix A. 

3.2 Module Construction 


Test panels were made using a new lot of Z-2591A polyurethane. 
This batch of the A part had a nonyellowing antioxidant. The 
panels were extremely slow in curing and showed bubble formation 
with time. After consulting with Development Associates/ a new 
lot of Part A was procured with a faster cure time. 

Additionally, the tanks containing the urethane components were 
plumbed into the house vacuum system. This storage under vacuum 
was successful in alleviating problems of bubbling. Several 
test coupons were made and showed no signs of bubble formation. 
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Two polyurethano modules were oonstructed using Tedlar as a top 
cover. The first module showed no signs of urethane-produced 
bubbling, however there were bubbles present which had come from 
beneath the cells. 

To avoid problems with bubbles trapped beneath colls, a modified 
process for assembly of the polyurethane layup wa>? developed. 

The procedure follows: 

1. The hardboard was prepared with wood ribs epoxied to the 1/8" 
4' K 4' board. 

2. Strips of 20 CP 3120 polyester film (3M) were glued to the 
back of the substrata using 4910 pressure sensitive adhesive 
(3M) . 

3 . The top of the substrate was spray painted with Krylon 
interior/exterior enamel 1502 (Borden) . This is an alkyd- 
based flat white paint. 

4c Holes were made in the substrate for electrical termination. 

5. The substrate was placed ribs down and a layer of .003 
Craueglas 230 placed upon it. 

6. The surface was flooded with polyurethane. 

7. The primed cell circuit was then placed on a board, one end 
of the circuit lightly held, and the board slowly pulled 
away to allow the circuit to fall on the flooded surface. 
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B. AtUiitional polyuroi haiu' waa thon pvuupod on the layup until 
tl\e surface was eonipletely eov#ered» 

‘1, Another layer of .00.1“ Oranevilaa was placed on the surface. 

10, rriwed (Dow Corninu *50020 in methanol) Tedlar lOOBG 
30UT was then rolled onto the layup to complete the modulo, 

11. A weiqhted steel plate was placed on the module and the. 
polyurethane allowed to eure. 


i'his nuHlule was the most success lul to date. Only a very few 
bubbles marred its appovirance. One problem was that bubbles 
became trapped under the first l*iyer of craneglas. This will 
bo overcome in the luture by first flooding the surface with 
urethane and then applyinu the Craneplas. This technique was 
very successful for the second layer. 


Section 4 . 0 


CONCLUSIONS AND RECOMMENDATIONS 


There are no conclusions and recommendations for this period. 


Section 5.0 
PLANNED ACTIVITIES 


During the next period construction of the qualification modules 
will continue. The thermal, structural, and optical final reports 
will be completed. 
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SUt!MAK\’ 


or rOOR QUArTTY 


the chrunt of the elettclcal Isolation Mnalyais petfomed in this phase 
of the encapaulatiop contract was to develop a luathod for avaluating the multi^ 
dimenaional affects of the module conductor geometry on the electric field* 
the electrical atreaa intenBification at edges and corners of the soler cells 
is not predicted hy the simple aeries cepacitance model developed in the first 
phase of the contract, the multidimensional electric field 'was calculated 
with the finlto-eleraent-hased NAStRAN thermal Analyeer, using the analogy 
between thermal and electrostatic fields* To evaiuete the maximum electric 
field accurately, the NASTRAN output had to he pout-proceaaed} a polynomial was 
fit to the predicted values of potential as a function of distance from the 
cell surface by the method of least squares, and the surface field was cal'« 
culatad from the heBt^fit polynomial, The accuracy of this method was verified 
by comparison with an exact solution for a geometry similar to that of typical 
solar cells. Two sample geometries were analyted — a square test coupon used 
in electrical isolation tests performed earlier in this contract and a fatilly 
of disc~shaped solar cells* these sample analyses demonstrsted tbst this 
finite-element method is a uaeful design tool for evaluating candidate module 
encapsulation designs* they elso showed geometric limits for which the model 
incurs numerical difficulties; 1) cells having very sharp edges, 2) cells 
much thinner than the dielectric pottant layer* and 3) cells much tUicKer thsn 
the dielectric pottant layer, these limits* some of which occur with 
advancad thin-film device modules* would require larger, specially designed 
computer models. It is recommended that electrical isolation analysis models 
suitable for advanced thin film device modules be developed, 
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Electrical safety is an important factor in the design of terrestrial 
photovoltaic modules, in which solar cells may operate at potential differences 
of several thousand volts relative to ground. Electrical breakdown could be 
a serious hazard to people (for example, to a worker cleaning a module) and to 
the hardware. Modules are required to be able to withstand at least 3000 Vdc 
without the occurrence of breakdown (References 1 and 2) . 

Short-circuit paths might be, for example, the pottant and organic front 
cover material in a substrate module, the pott^ut and wood-product /load-bearing 
member of a wood substrate module, or the pottant between the cell and metal- 
foil back cover of a glass superstrate module. In addition, manufacturing 
defects such as cracks and bubbles will have a deleterious effect on the 
ability of a module to withstand high voltage stresses. Therefore electrical 
isolation is a key factor in the design of the module encapsulation system. 

Accordingly electrical isolation has been investigated throughout this 
encapsulation contract. In the first part of the contract, a simple series 
capacitance model (shown in Figure 1) and a simple series resistance model 
were developed for determining the electric field strength in each material 
layer of the encapsulation system (Reference 3). These calculations were 
employed to specify the minimum layer thicknesses required to keep the electric 
field in each layer lower than the dielectric strength of the material in 
that layer (assuming that there are no flaws or manufacturing defects). In 
the second part of the contract, the breakdown voltage was measured for 
specially designed test coupons employing simulated solar cells (see Figure 2) 
and a variety of encapsulation schemes, and these measurements were compared 
with the predictions of the series capacitance model (Reference 4) . The key 
results were as follows (see Table 1) : 

r •- , 

i - , * 

s For most of the encapsulation schemes tested, the breakdown 
voltages predicted for a particular sch^e were within the 
range of the values measured for the test coupons of that 
scheme. 


1.0 INTRODUCTION 


1.1 BACKGROUND 
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ANALYTICAL 
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DESION CniTEniA 

• E| < Sj rOR ALL LAYERS 

E| * ELECTRIC FIELD STRENGTH IN 
LAYER 1 

S| - DIELECTRIC STRENGTH 

• Vp - 3000 VOLTS DC (LSA REQMT) 


Figure 1. Series capacitance model developed in first part of contract for determination of 
electric isolation characteristics of photovoltaic modules. 



Figure 2. Test set>up for electrical isolation tests performed in 
second part of contract; dimensions in inches. 
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TABLK 1. COMPARISON Of moiCTED AND MEASURED 
VALUES OR At BREAKDOWN 





'■o 

at Breakdown, 

kV 

Coupon 

type 



Measured 

... 

Predicted 
with Series 
Capacitance 
Model 



Sire 

Description 

Max. 

Min. 

Mex. 

Min. 

A 

Rront 

A mil tediar, 18 mil EVA/GG* 

19 

10 

U 

ii .8 


Back 

18 mil EVA/CG, 1 mil 
A1 Polyastere 

11 

1 

11.2 

e .e. 

B 

Rront 

1 mil Tediar, 18 mil EVA/CG 
36 rail EVA/CG, 1 mil 
A1 Folyaster 

19 

12 

11.4 

11,3 


Back 

13 

5 

22.3 

n 

e^«, « 

jp% 

V. 

Rront 

1 rail tediar, 18 mil EVA 

21 

5 

11.4 

11.3 


Back 

18 sil EVA/CG, 125 mil wood 

25 

8 

24,3 

23.9 

Dl 

Front 

1 mil tediar, 36 mil EVA/CG 

25 

12 

22.6 

22.5 


Back 

36 mil EVA/CG, 125 mil wood 

25 

21 

26.8 

26.0 


Front 

1 mil tediar, 36 mil EVA 

23 

10 

22.6 

22,5 

m 

Back 

36 mil EVA/CG, 125 mil wood 

25 


26.8 

26.0 


*BVA/CG • EVA with Ccai>egia» 

A1 Vdiyeafeer - Aiuwinisah Poiyascer 


• On cl»tt ochw hAn4, chn breakdowu^voitage rangea much 

wlday Chan thoaa ptadiccad by aubacitution o£ the known mateciai»> 
property rtngei into the aeriea capacitance model. 

• ttw measured breakdenm volfeagea had a weaker dependence on the 
potcant thickness than predicted by the series capacitance 
model. 

w All the electrical Allures occurred at the edges of the simulated 
solar cells in the test specimens. 

these tesulta indicated the need lor a wore realistic model than the 
series capacitance model, the series capacitance model assumes chat the elec- 
tric field is unifor»\ ii^ a layer , thac is, that the lines of force are straight 
















and Cbac the aqulpotontiala are equally apaced parallel planes. Actually the 
equlpotentlala crowd together at aliarp edges and cornera of the cell, as ahown 
In Figure 3. Gonaequently the electric field, which is the gradient of the 
potential, has its maximum magnitude at auch locations. In Che absence of 
flaws, the onset of breakdown can be expected to be governed by and to occur 
at auch sharp edges and corners. This conclusion is well known (Reference 5) 
and is consistent with the results of the aforementioned tests performed in 
this contract. Therefore a realistic model must account for the multidimen“ 
sional effects of the module conductor geometry on the electric field. 



Figure 3. Electric field distribution in module. 
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1.2 UEAT TRANSFER ANALOGY 

Fortunately three-dimensional electrostatic fields can be calculated by 
available heat transfer computer programs, because of the analogy between 
thermal and electrostatic fields (see Figure 4). Laplace’s equation governs 
both the temperature distribution In conduction heat transfer and the poten- 
tial field In electrostatics. Tlie electric field, the quantity of Interest In 
the electrical isolation analysis. Is analogous to the temperature gradient. 
Because of this analogy, L. B. Duncan of Hughes suggested using the NASTRAN 
Thermal Analyzer computer program to calculate the electrlc-f leld enhancement 
resulting from candidate encapsulation schemes. 


• DIFFUSION EQUATIONS 

r«*rV+p -0 EUECTBICAL 

V.ttTT>4-pc(|I) thEHMAU 


THEl^MAL FASAMETES 

ELECTRICAL FARAMETER 

TEMFEFIATUSE 

■■ 

FOTENTIAL 

V 

THERMAL CONDUCTIVITY 


FERMITTIVITY 

f 

temferature gradient 


ELECTRIC FIELD 

T- “TV 

HEAT FLUX 


ELECTRIC OI8FLACEMENT 

"o-tT 

INTERNAL HEAT GENERATION 

■ 

CHARGE DENSITY 

p D 


Figure 4. Thermal-electrical anilogv. 


1.3 PRESENT INVESTIGATION 

Accordingly, an investigation of the feasibility of performing electric 
field analyses for photovoltaic module geometries was performed. This inves- 
tigation is the subject of this report. 
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The specific tasks were as follows; ^ POOR QUALfTY 

1) Review the literature on previous electrical analyses performed 
with finite-difference or finite-element computer programs using 
the heat transfer analogy. 

2) Formulate and solve multidimensional sample problems to check 
the aforementioned series capacitance model and to investigate 
edge and corner effects. 

The literature search showed that the finite-element method, employed in 
NASTRAN to solve the partial differential equations governing thermal and 
structural phenomena, has been employed to calculate electric fields since 
the late 1960s. Books have been written on the method (for example. Refer- 
ence 6) , and the MSC/NASTRAK program used for thermal and structural analyses 
at Hughes has been employed for electrical analyses (References 7 and 8) . 

These findings Indicated the likely success of the application of this method 
to the photovoltaic module electrical Isolation analysis. 

The remainder of this report describes this application. Section 2 
describes the development, verification, and use of the method. Section 3 
describes the sample geometries analyzed and the results. And Sections 4 and 
5 present conclusions and rocommendations. 

2.0 Mbfb'OD 

2.1 NASTRAN THERMAL ANALYZER 

The finite-element-based NASTRAN Thermal Analyzer (Reference 9) is a 
general-purpose heat transfer computer program. As an integrated part in the 
NASTRAN (I^A structural analysis computer program) system, this thermal 
analysis capability is fully capable of rendering temperature solutions and 
heat flows in solids subject co various boundary conditions, which range from 
prescribed temperatures at grid points and specified thermal loads to convec- 
tive and radiative modes of heat transfer at boundary surfaces in both steady- 
state and transient cases. This heat transfer computer program has been 
developed by an application of existing functional modules in NASTRAN, which 
were designed originally for the purpose of structural analysis, and the addi- 
tion of new modules including new elements and new solution algorithms. 
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2.2 CALCULATION OF MAXIMUM fIeLD .Of fOORiQU" 

w» • 

The quantity of Interest for the photovoltaic nodule electrical isolation 
analysis is the maximum electric field in the module. For the geometry shown 
in Figure 3, the maximum field occurs at the cell edge or corner, and the field 
vector there is perpendicular to the cell surface. For a finite edge or 
corner radius of curvature, the precise location of the maxlTBUn field is not 
known a priori. The problems attacked here were to locate and evaluate the 
maximum field from the NASTRAN output, namely, the average potential In each 
of the analysis elements. 

This endeavor was not straightforward. NASTRAN calculates the field 
simply by 


^ASTRAN " » 

where AV ■ pstenclal difference between adjacent analysis elements and 
d ■ spacing between these elements. 

This formula significantly underestimates the surface field because, as seen 
in Figure 3, the potential if, a nonlinear function of the distance D from the 
cell surface. In otiter words, the function V(D) has a non-zero second 
derivative. Trial computer runs indicated Chat the field computed by NASTRAN 
underestimated the surface field and that the use of a finer computational mesh 
or modifications in the shape of the mesh did not solve the problem. 

The successful solution was to: 

1) Calculate the potential distribution with a moderate size 
(around 300 elements) NASTRAN model, 

2) Perform a least-squares curve fit to the NASTRAN output values 
of the function V(D), and then 

3) Calculate Che surface field from the best-fit formula. 

The development of this procedure Is described In Appendix A. 

2.3 VERIFICATION OF METHOD 

A geometry, similar to those of interest for photovoltaic applications, 
for which there is an exact solution of Laplace’s equation was discovered. 
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This g«oio6tr>‘ coneistB of i«\ infinitB«iin«lly thin pl«n«r electrode (V - 1) 

located halfway ^ftweec twQ ground planes. The hlgh**volcage electrode extends 

to infinity toward the right* and the ground planes extend to infinity toward 

* » 

the left and the righ^. The equipotential surfaces ssynpcotically approach 
parallel planes and can be destiribed as blunted knife-edged slabs. The 
details of this exact soluation are described in Appendix B. 

The equipotentials reseable the sliapes of solar cells and thus can serve 
as A model of solar cells tvsving axial symmetry (discs) or cells having planar 
geometry (strips). Any equipotential can be thought of as the surface of a 
solar cell. 

To verify the aforementioned method for calculating the maximum electric 

field* a NASTRAN model of the blunted knife-edged slab was developed for the 

region between the V •• 0.6 equipotential aurfnee and the ground plane. This 

aurfacc and two other equipotentials are plotted in Figure 5. The geometrj' 

corresponds to a ratio of cell thickness to pottant thickness having the value 

t /t ■ 0.5. The KASTRAK model had 26A elements and 308 grid points. The 
^ P 

maximum field occurs at the surface halfway between the ground planes. 

The maximum fields predicted by the exact solution and by the best-fit 
polynomial to the HaSTRAN output agreed to within 3 decimal places, which is remark- 
ahly close agreement. The potential distributions predicted by the exact solu- 
tion and by NASTRAN also agreed very closely* as shown in Figure 6. 

This excellent agreement with the exact solution verifies the accuracy of 
the method described in the previous section for calculating the maximum field. 

It also indicates that a moderate-sire model consisting of a few hundred analysis 
elements provides sufficient accuracy for typical geometries of interest. 

2.4 USE OF METUOP 

This method can be used in providing for the electrical isolation of photo- 
voltaic modules as follows. To avoid breakdown, the module should be designed 
so that 

E < S, (1) 

m ’ 
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Figure 5. Blunted knife-edged $l»b, exact solution (tj,/tp - 0.5). 
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POTENTIAL MIDWAY BETWEEN GROUND PLANES, V (u . 


where 


in 
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• • h « 

maximum electric field in module and 


S dielectric strength of pottant. 


The maximum field is a function of the following parameters: 

• The voltage difference Vq between the solar cell and the 
ground plane (the modules are required to withstand a voltage 
difference of at least 3000 Vdc without breaking down), 

• The thickness t^ of the pottant, and 

• The geometry of the module (specifically, the effects of 
sharp edges or corners on ccnductors) . 

In the previous series capacitance model, the electric field in each 
layer was assumed to be given by the value for infinite parallel plates: 


PP 


V /t . 

0 P 


( 2 ) 


(This is a good approximation to the field except near edges and corners of 
conducting surfaces •) This enabled the minimum allowable pottant thickness to 
be calculated by 


E < S, (series capacitance model) (3) 

PP 


which yields 


tp > V^/S. (series capacitance model) 


(4) 


With the present method, we define y as the ratio of the field for 
infinite parallel plates to the maximum field: 


E V /t 

PP , o P 

E E 

m m 


(5) 
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The quantity y, which is a function only of the module geometry, is an 
electrical stress intensification factor. It indicates the amount by which 
the module's voltage^withs tending capability is reduced because of sharp 
edges or comers. For example, If y “ 0.7, the maximum field in the module 
would exceed the pottant's dielectric strength at a voltage difference 70% 
as large as would be predicted by the series capacitance model. Thus the method 
described herein for calculating the maximum field can be employed to evaluate 
the electrical isolation characteristics of candidate module geometries, as well 
as to specify the minimum allowable pottant thickness. The design criterion, 
replacing (4), is 

y t > V/S, (6) 

p o 

This method does not predict the voltage at which breakdown will occur 
across the entire pottant layer for a given module. Rather it predicts the 
voltage at which breakdown could begin in a microscopic region. Predlcticn 
of the development of this onset of breakdown into a conductive path across 
the entire pottant layer involves complicated phenomena (see Reference 5) and 
would be a big step beyond the present analysis. The work reported herein 
was limited to the method for predicting the electrical stress intensification 
factor y. 

3.0 ANAX.YSIS OF SAMPLE GEOMETRIES 

3.1 SQUARE ELECTRICAL ISOLATION TEST COUPON 

The first sample geometry analyzed was a typical test coupon used in 

the electrical isolation tests performed earlier in this contract and reported 

in Reference 4. The back side of coupon Type A (see Table 1) was modelled. 

The copper test coupon and the finite-element model of the comer region are 

shown in Figure 7. The analysis model consisted of 520 elements and was 
constructed so that the radius of curvature of the corner in the plane of the 
coupon could be varied as a parameter,, The edges of the coupons were assumed 
to be infinitely blunt. (The radii of curvature of the corners and edges of the 
actual test coupons were never measured, and no information on typical values 
was available.) 


- 25 - 



OTTANT 




tiWRUNAic pm: t:/ 

Of poon QUAirrv 

Th« d^ecbrlcnl IntisnRiflcAtlon fndtot y 1 h plotted veteuii 

the ratio of the radlua of curvature of the teat coupon to the pottant thlck*^ 
neaa in figure 8. The cur\*e appeiira to behave correctly in the liroltas 

• y 0 »a 

tl\ac 1 b, the electric field becoinea infinite in the liodt of a perfectly aharp 
corner 

• y 1 «B R/t **• , 

p 

tliat in, the electric field approaches the value for parallel platen in the 
limit of a very blunt corner, 



p|0ure 8< Spamv «*« coiipijn »» eifnGl oi coupon cornier rfcdiui of curv»tur« 
on flldcuicttl itnsii tntfiniHlcJUlon fi«cior. 
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The test data listed In Table 1 are for the voltage at which the entire 
layer broke down, not the voltage at which the breakdown began In a saiall 
region, and are not comparable with the predictions of the present analysis. 

3.2 DISC-SHAPED SOLAR CELL 

The other sample geometry analyzed was the family of disc-shaped solar 

cells shown In Figure 9. The radius of curvature R of the sun-side edge was 

assumed to be the same as that of the other edge. The electrical stress 

Intensification factor y Is a function of the parameters R/t^ and x(*t /t^), 

c p c 

where t « cell thickness. As shown In Figure 9, R/t varies between the 
c c _ 

values 0 (the limit of a perfectly sharp edge) and 0.5 (the limit of a 
semicircular edge having no flat portion). The scale drawings below the 
abscissa In the figure Illustrate the edge shapes for various values of R/t^. 

This family was selected for analysis because It Is a good model .of 
present and future solar cells. Of particular Interest were the following 
limits: 

• R/t -► 0. 

c 

This Is a sharp-edged cell (see Figure 10) and represents the limit of the 
type that Is produced currently by single-crystal growth and laser scribing. 
On physical grounds, we expected that y 0 as R/t^ 0, that Is, the maximum 
electric field Is Infli^te for a perfectly sharp edge. 

• t /t -*• 0 for R/t • 0,5. 

c p c 

This Is a thin, blunt-edged cell (see Figure 10) and represents advanced 
thin-fllm devices. We held R/t^ -0.5 because it seems Implausible that the 
edge of a thln-fllm cell will have a flat portion. 

• X (»t /t ) -► 0 for R/t ■ 0.5. 

pc c 

This Is the limit of a very thick, blunt cell or the limit of a veTry thin 
dielectric layer. We expected on physical grounds that y -*• 1 as x 0, that 
Is, the maximum electric field approaches the parallel-plate value for an 
infinitely blunt cell. 
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FI()Ur* 10, Typioil plots of equipotantials for ditc-shapad cell, 


The analysis model consisted of approximately 400 elements and 460 grid 
points. It was constructed so that R/t^ and x could be varied as parameters. 

Figure 9 shows the predicted dependence of y on R/t^ for several values of 
X. In accord with our aforementioned expectation, the computed curves (shown 
as solid curves) have been extrapolated as dashed curves to the origin. As 
explained In Appendix A, the model Incurs numerical difficulties for small 
values of R/t^ (sharp edges) and large values of x (thin cells). Some of the 
results In these regimes appeared physically wrong and were not used. The 
model could be extended to these regimes In the future. 

Figure 11 shows the dependence of y on for R/t^, “ 0.5 (a semicir- 

cular edge) . The scale drawings below the abscissa show the geometry for 
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Figure 1 1 , Disc*shaped cell ^ effect of cell thickness on electrical stress intensification factor. 


various values of t /t . For the thinnest cell analyzed (t^/t « 0.02 (x ■ 

30)), y ■ 0.19; that is, the maxltnum electric field is approximately five times 

as high as for parallel plates having the same spacing. For such a small value 

of V /t , the model has numerical difficulties, and it is difficult to deter- 
c p’ 

mine the limit of y as t /t -^0. This limit, which is Important for advanced 

c p 

thin-film solar cells, could be evaluated by further work on the model. 

To determine the limit as x •+• 0, we plotted y versus x in Figure 12. As 

was mentioned previously, we expected on physical grounds Chat y 1 as x 0. 

Calculated values for x ^ 0.5 (t /t £ 2) do indeed appear to lie on a curve 

c p 

that is approaching the (0,1) point. Me also computed the value for x 0.25 



4), and it falls near this curve. However, the present model begins 
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Figure 12. Diic-«Hsped cell — llmilt of electrical stress intensification 
factor as pottant becomes very thin. 
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to incur some numerical difficulties at x « 0.25, and it would take some labor 
to develop a model for this regime. We do not think that this is warranted, 
because accurats values of y for large values of t^/t^ can be obtained by extra- 
polating the curve to (0,1). To demonstrate this, we fit the computed values 
for 0.5 X ^ 20 to a 5th degree polynomial. (We ignored the values for 
X > 0.25 and 50, for which the model is less accurate.) The best-fit formula 
was 

y - 0.99279 - 0.409946 x + 0.140055 x^ 

- 2.46179 X 10"^ x3 + 1.87591 x lO"^ x^ 

- 4.74338 X 10“^ . (7) 


The index of determination was essentially unity, indicating that Eq. (7) 
is a very good fit. The curve resulting from this formula is smooth and 
monotonic for values of x as large as approximately 2. Therefore it can be 
used for predicting values of y not calculated by the finite-element model for 
X _< 2. For values of x much larger than 2, the best-fit curve is oscillatory, 
characteristic of a high-degree polynomial. The important result for the 
present purpose is that the best-fit curve passes very close to (0,1) (actually 
it passes through (0,0.993)), even though we did not constrain the best-fit 
starve to pass through (0,1). This is strong evidence to support our hypothesis 
that the actual curve passes through (0,1). 

The quantity of Interest for the isolation design of a module is the 
maximum electric field E , which is given by Eq. (5). It can be written as 


E 


m 



( 8 ) 


The mavHniiwn field can be evaluated from the finite-element results as follows. 
For a given applied voltage and cell thickness, E^ is a function only of the 
product yx, which is a function only of the module geometry. As an 
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example, consider the family of semlclrcular-edged solar cells shown In 
Figure 12. For this family, y Is a function only of x. The product 



( 9 ) 


Is plotted versus x In Figure 13. It Is seen to be a monotonlcally Increasing 
function of x for x ^ 20. Therefore the maximum electric field decreases with 
Increasing pottant thickness, for poctant layers as much as 20 times thicker 
than the cell. It Is recommended that this effort be continued to extend the 
finite-element model for values of x ^ 20, which are of Interest for advanced 
thln-fllffl devices. For comparison, the predictions of the series capacitance 
model also are shown In Figure 13. Because of the electrical stress Inten- 
sification at the cell edges, the pottant has to be much thicker to reduce the 
maximum field to a given value than Is predicted by the series capacitance 
model. 


4.0 CONCLUSIONS 

1) A finite-element method for evaluating the electrical-isolation 
characteristics of photovoltaic modules has been developed. 

2) Its accuracy has been verified by comparison with an exact 
solution for a geometry similar to that of typical solar cells. 

3) The utility of the finite-element method as a design tool has 
been demonstrated by the analysis of two sample geometries. 

4) Accurate calculation of the maximum electric field required careful 
analysis. The finite-element models had to avoid large elements 
adjacent to small ones and had to avoid elements with large 
aspect ratios. The NASTRAN output had to be curve fit to 
calculate the maximum field accurately. 

5) Larger, specially designed computer models would be required for 

the following three limits: cells having very sharp edges, cells much 
thinner than the dielectric pottant layer, and cells much thicker 
than the dielectric pottant layer. 


5.0 RECOMMENDATION 

The capability of the finite-element model should be extended to the 
geometries of advanced thln-fllm device modules. 
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Figure 13 . Diic-thaped cell - effect of pottent thicknew on mexlmum electric field. 
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APPENDIX A 

DEVELOPMENT OF NUMERICAL METHOD FOR CALCULATING MAXIMUM 

ELECTRIC FIELD 

This appendix elaborates on the description in Section 2.2 of the method 
for calculating the maximum electric field. 

Figure A-1 ahovs the output of a typical NASTRAN run. It consists of the 
potential at each grid point In the finite-element model. For this geometry, 
the maximum potential gradient occurs on the surface normal denoted by grid 
points 1, 2, 3, and U. The maximum field, which occurs at the cell surface 
(grid point 1), is given by 


E 


m 


dV(D) 


dD 


D-0 


(Al) 


The actual function V(D) (the solid curve in the figure) is seen to have a non- 
zero second derivative. NASTRAN, however, calculates the field simply by 
fitting a straight line through adjacent grid points: 


m 


”2-»i 


(lUSTIUll) 


(A2) 


Comparison of the two dashed lines in the figure showed that the straight- 
line approximation underestimated the surface gradient by 61. The difference 
between the two values of the surface gradient was 10-201 for many of the other 
geometries analyzed in this project. The error in the surface gradient result- 
ing from the straight-line approximation was difficult to detect by visual 
Inspection of the graph, but it shoved up when surface gradients were cal- 
culated. To improve the accuracy, we tried 1) using a finer computational 
mesh and 2) modifying the shape of the mesh. Neither modification yielded 
adequate accuracy, because Eq. (A2) inherently underestimates the surface 
gradient . 
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Figure A-1. Underestimation of maximum field by NASTRAN 
illustrating need to post*process output data. 
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Accurate values finally were obtained by poat-proceaslng the NASTRAN 
'output. The first step was to determine the direction of the maximum surface 
gradient by Inspection of the HASTRAK output. Then the Ath-degree polynomial 
that was the best fit to the function V(D) along this direction was determined 
by the method of least squares. The curve fit used the potentials at four or 
five grid points, depending on the geometry. Then the surface derivative was 
evaluated from the best-fit polynomial and equated to the maximum field, per 
£q. (Al). With this post-processing, a moderate-size (around 500 elements) 
NASTRAM model provided sufficient accuracy, except for extreme geometries 
discussed in the next paragraph. As mentioned In Section 2.3, the maximum 
field predicted by this method agreed very closely with that calculated from 
the exact solution for a geometry for which an exact solution of Laplace's 
equation exists. 

Even with this method, however, care Is required In developing the finite - 
element model to obtain accurate results. Analysts should pay special attention 
to 1) the fineness of the computational mesh and 2) the shapes of the finite 
elements. A fine mesh Is required In regions having large gradients, such 
as near edges and comers. But the analyst must take care to gradually make 
the mesh coarser with Increasing distance from the edge or comer. We found 
that a model having a large element adjacent to a small one gave Inaccurate 
results. MASTRAK also Incurs numerical difficulties with models having 
elements with aspect ratios greater than about 1,500:1. Therefore elements 
with large aspect ratios should be avoided. For these reasons, models of 
1) edges or comers with very small radii of curvature, 2) very thin cells, 
or 3) very thin dielectric pottant layers must be very large and thus 
expensive. 
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APPENDIX B 

EXACT SOLUTION FOR BLUNTED KNIFE-EDGED SLAB 


The geometry for which the exict solution, described In Section 2.3, exists 
Is shown In Figure B-1. The geometry Is planar and Is specified by Cartesian 
coordinates u and V. The ground planes are located at v - 0 and v • ti. The 
other electrode, at which V ■ 1, Is located at v “ ir/2 for u ^ 0. 


t 

ir 


X'2 


V - 1 


V - 0 


0 


Figure . Configuration for exact solution, 


The potential function In the region between the ground planes can be 
determined by conformal mapping, using the Schwarz-Chrlstoffel transformation. 
The solution Is (Reference Bl) 


T 



2 1/2 
8 ) 


/s, 


where 


T ■ tan nV, 


(Bl) 


(B2) 
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B - - 1 + (1 + 2c”^“ 


C 4- e 



1/2 


» 


(B3) 


C ■ COB 2v, (BA) 

The quantity s is positive for V < 1/2, zero for V « 1/2, and negative 
for V > 1/2. 

The equipotential surfaces have the shapes shown in Figure B-2. They 

are blunted knife-edged slabs. The edge is located at the point (u , tt/2), 

o 

where 



2(1 4- 




[(1 4- t 2)^/2 4- sign(s)3. 


(B5) 


V 

♦ 10 , «r) 


I 

I 

I 



< 0 , 0 ) 


Figure B>2, Equipote^tials, 


Each surface approaches a pair of parallel planes as u ®, that is, 

lim T “ 4 tan 2v. (B6) 

u-+®> 
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11m V ■ S 


u 


It au - Tt 

u-'jV, V>.J 


V < 


(B7) 


Thus an equlpotential surface can be thought of as a solar cell having a 
thickness 


^c “ 


(,-fv) 


V - tt(1 - V) 


<B8) 


encapsulated in a dielectric pottant having a thickness 




(119) 


Therefore 


t v 

p 


(BIO) 


The value V “ 1/2 corresponds to t /t ■ 2. The electric field approaches the 

c p 

asymptotic value 




11m E (upv) • 

V Cn tl 

U-H» P 


(Bll) 
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The nagnitude of the electric field vector le given by 


2 2 
E » {(3V/Su)^ + OV/3v)‘^l 


(B12) 


By evaluating these derivatives, the electric field can be written simply as 


£ 


2/r 


(s + V) 


1/2 


(B13) 


The maximum field E is determined as follows. For V < 1/2, s > 0, 

m — — 

Therefore the maximum field occurs where s ■ 0. From Eq. (B3) , s • 0 at 
u ■ “>. Therefore E^ • E^ and is given by Eq. (Bll) . In other words, the elec- 
tric field does not exceed the parallel-plate value in the region where V ^ 1/2. 
For V > 1/2, s <. 0 * Therefore the maximum field occurs where $ has its 
minimum value. From Eqs. (Bl) - (B3) , we have 

s - 211 + (1 + T^)^^^C]/T^. (B14) 


The quantity T « tan rV is constant on an equipotential surface. The quantity 
C • cos 2v can vary between the values -1 and 1. The min im u m value of s 
occurs where C » -1 (that is, at v ■ r/2) and is given by 


s - 2[1 - (1 + t2)^^^]/T^. 
m 


Therefore the maximum field is given by 

2/tt 




2/ti 




1/2 


. V < 1/2 
, V > 1/2 


(B15) 


(B16) 
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The electrical stress Intensification factor y, defined by 


y • ’ 


Is given from Eqs. (Bll) and (B16) by 


, V < 1/2 


(8^ 1)^^^ , V > 1/2 

in 


It Is plotted versus in Figure B-3. 
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Figure B*3. Electrical stress intensification fector for blunted knife-edged slab {exact solution). 
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